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Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds 
have been developed recently. For the proposed flight conditions the air in the shock layer that 
envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational ex- 
citation, and ionization. However, these processes occur at finite rates which, when coupled with 
large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non- 
ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. 

A numerical method to solve the equations that describe these types of flows in two dimen- 
sions has been developed. The state of the gas is represented with seven chemical species, a sepa- 
rate vibrational temperature for each diatomic species, an electron translational temperature, and a 
mass-averaged translational-rotational temperature for the heavy-particles. The equations for this 
gas model are solved numerically in a fully coupled fashion using an implicit finite volume time- 
marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and 
flux-dependent differencing is employed to maintain stability. 

The numerical method has been tested against several experiments. The calculated bow shock 
wave detachment on a sphere and two cones was compared to those measured in ground testing 
facilities. The computed peak electron number density on a sphere-cone was compared to that 
measured in a flight test. In each case the results from the numerical method were in excellent 
agreement with experiment. The technique has been used to predict the aerothermal loads on an 
Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the 
current physical model of high temperature air is appropriate and that the numerical algorithm is 
capable of treating this class of flows. 
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Research Objectives 


r» i-j 

c~h a 
IK B 


D o3 

•5S O 


o o 


502 


Demonstrate feasibility of numerical solutions to reacting flows 
o stiff equations 
o large source terms 



Justifications 
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effects of thermo-chemical nonequilibrium 
develop new modeling techniques 



Technical Approach - Thermodynamic Model 




Technical Approach - Conservation Equations 
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Technical Approach - Equations of State 
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Technical Approach - Chemical Model 
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use two-temperature reaction rate (TT V ) model of Park 



Technical Approach - Numerical Algorithm 



stability in subsonic and supersonic regions 
capture bow shock wave 



Validation - Peak Electon Number Density 
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Results - Cylinder 


• Interferogram of 2 inch diameter cylinder in N 2 and N . 
tioo = 5 .59 km /s , Moo = 6 . 1 , Re = 12000 ,-0 = 5.5. 
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OF POOR QUALITY 
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Status of Computer Code 



other gas models 



Conclusions 
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